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INTRODUCTION 


The  extormil  design  of  most  projectiles  was  established  during  and  prior  to  World  War 
11.  hew  configuration  changes  aimed  at  increasing  the  range  of  these  projectiles  have  been 
incorporated  since  that  time.  To  achieve  maximum  range,  the  total  drag,  which  includes 
pressure  drag,  skin  friction  drag,  and  base  drag  must  be  minimized.  In  1969  one  of  the 
authors*  developed  an  improved  projectile  shape  by  deriving  an  approximate  minimum 
total-drag  projectile  based  on  a large  amount  of  extu'rimental  data.  Although  a number 
of  theoretical  analyses  has  determined  the  shape  of  minimum  drag  bodies  by  considering  the 
pressure  drag  and  in  .some  cases  skin  friction  drag  also,  none  has  included  the  base  drag  as 
well.  This  report  develops  an  analytical  method  for  theoretically  predicting  the  projectile 
shape  which  yields  the  minimum  total  drag  for  a fixed  length,  diameter,  and  supersonic 
Mach  number  ( 2 ^ ^ 6). 


A review  of  the  literature  revealed  extensive  research  on  minimum  or  low  drag  shapes, 
but  only  part  of  the  total  drag  was  minimized.  References  2 through  10  all  deal  with 
bodies  of  low  or  minimum  wave  drag  (pressure  drag).  Of  particular  interest  is  the  body  of 
minimum  wave  drag  developed  by  von  Karmatt.'  The  presstire  distribution  was  estimated 
from  slender  body  thei>r> . and  then  the  wave-drug  integral  was  minimized  by  the  calculus 
of  variations,  't  he  resulting  Innly  sl»a|H'  is  given  by  the  relation 


r a 


] sin  2t^ 


where  ^ ^ cos  * ( 1 2vA  ).  Stivers  and  Spencer*  * eulculaied  the  total  minmunnHirag 
characteristics  of  IVnir  fa imlies  of  slender  Innlies  at  Mach  numbers  from  2 to  12.  Of  these 
four  families,  they  found  that  the  Searvltaack  pridile  |»rovides  the  lowest  total-drag 
coefficients  at  zero  incidence.  I hey  did  not.  however,  eoiisider  prid'iles  outside  these  four 
families.  Iletisc,  there  is  no  reastm  to  IvUeve  that  the  Searvllauck  Innly  is  the  optimum 
sltajv. 


( ole*'  alw»  attempted  trr  tiiid  a Nuly  »rt  ittmiinutii  wave  drag  by  using  the  Newtonian- 
UuH‘tiiann  theory  lor  slender  IwHhes  to  estimate  the  prevsutv  distribution,  llte  body  that 
tils  these  conditions  is  the  2t%(-|Hnver  law  iHHly  which  rs  descritred  by 

I lie  von  Karman  binfy  has  a |Htinied  nose  wliervas  ibe  2/.Utwiwvr  law  IhhIv  has  a small 
blunk'd  nose.  In  .idsitium  to  iK'gleciing  skm  irieiion  drag  and  Inise  drag,  IhmIi  i»f  these 
sliaires  wvre  dvwvto(vd  using  slender  iHMiy  ilu’otisw 


Mielc**  attempted  to  optimize  the  body  sliape  by  minimizing  the  sum  of  the  pressure 
drag  and  skin  friction  drag  for  slender  bodies.  As  in  previous  cases,  the  body  shape  is  simply 
an  ogive  since  base  drag  was  not  considered.  However,  in  the  Mach  number  range  at  whicli 
many  present  and  future  projectiles  travel  (0.7  < < 5.0),  base  drag  is  quite  important 

and  represents  a significant  portion  of  the  total  drag.  In  the  present  report,  approximate 
yet  reasonably  accurate  techniques  are  used  to  predict  pressure,  skin  friction,  and  base  drag. 
'Mien  the  calculus  of  variations  is  applied  to  determine  the  body  shape  whiclr  gives  minimum 
total  drag  when  the  length,  diameter,  and  Macii  number  are  held  constant. 


ANALYSIS 


The  Mach  number  region  of  interest  here  will  be  primarily  1 ^ ^ u 1 his  covers 

the  range  of  interest  for  projectiles  with  ininiinum  time  of  llight  rei|uirements  such  .ii> 
those  used  for  antiaircraft  applications.  It  also  overlaps  the  Maeh  number  range  ot  interest 
for  most  conventional  projev tiles  in  use  today.  0.5  r*,  d/,  .V  lienee,  some  coinp.mson 
can  be  made  k'lween  the  theoretical  shape  derived  herein  with  the  semiempnical  slnipe 
derived  in  Reference  1 for  Mach  numbei\in  the  1 nv  supersonic  regime.  .Also,  since  it  has 
been  shown  that  the  Innly  of  ininimutu  drag  fur  fi.sed  length  and  diameiei'  Im.s  a •.lightly 
blimteU  nose  (bluntncss  ratio  **  0. 1 ),  only  blunt-nosed  eoniiguratioiis  will  be 

considered, 


The  approavii  here  is  to  use  appro.viinale  analy  iieal  teeliniques  to  calculate  the  wave, 
skin  friction,  and  base  drag,  Approximate  teclmuiues  aiv  necessary  since  theoptimi/ation 
procedure  involves  several  iteralrons  of  body  getmietry,  and  use  m mure  “exact”  numerical 
tecltnttiiies  would  make  tire  ctmiputaiioiial  lime  and  storage  prohibitive.  On  the  other 
hand,  approximate  leclinniues  wijl  K*  compared  with  im«e  exact  methods  to  determine 
their  accuracy, 


URACi 


The  lulul  drag  coctfivieut  for  a luHly  of  twoluiioit  at  zero  angle  of  attack  is 


S jf  r|/t v»|^  tt'v  * ^ 


tl) 


Uefetting  l*»  iTgure  1.  ( „tvi  is  the  prevMire  ctrenkieni  along  the  surface.  is  the  mean 
skm  friction  eoeflieleni.  and  is  the  base  pressure  e«KdTielenl  (C*,^  is  generally  negative 
whieli  makes  the  third  term  in  lw)uuiion  1 1 f positive),  t he  iifsi  term  on  the  right  side 
of  I quaiton  (1 ) is  the  wave  oi  pressure  dratt.  the  second  term  is  the  skin  IrietUm  drag,  and 


du  dp 


I it'urv'  I t IUhIv 


!lu'  tiuul  term  is  the  Imh*  dra^.  I liv  priilMetu  rs  it»  deUntnitK'  the  IhkIj  sha|v  vvhidi  nu«imi/vs 
I ijualiuti  < 1 1 suhieet  llte  sH*nsiramt\i»i  ti\eU  Inniy  knjMit.  il»amele»v.md  Maeh  mtmlvr, 

io  mmimi/e  the  ti»i.it  dr.it:.  metluHls  ate  twvded  t**  |>ie*l»wi  ( m . amt  ( as 
luitdruits  ul  IhhI>  slutv  ami  Mas li  mmdvr  liiv’sc  hh'UuhIs are  disvusH’d  ktterty  lK*Um. 


I'uretHHly  umJ  ?\i'U*rl>ml>  1‘ressure  l'»K*lTicu*i« 

I tis‘  rmelnnly  jirvssuts*  c*x*ntcu*»rl  is  caUulated  iK>m  itx*  iruidiitvxt  \s*\vhurum  j'livssurc 
disuilnitUm’ 


'dire* 


svlieis*  «*  IS  Uw  KhIv  slujv 


r'lvi  laiH* 


ami  the  stagnation  pa'ssiire  coelTidcnt  bcliind  a normal  shock  wave  is 


^’f)  " 71/-; 


(7  + I )A/i 


> 

1 T 


7 + I 


(7  I ) 


'} 


(4) 


i (liiation  ( 2)  is  used  to  calculate  the  pressure  eoerneient  iVom  the  hUinted  nose  up  to  tlie 
point  ormavinuiin  tliickness  where  d - 0.  Although  liiuation  1 2)  gives  C^,  = Oat  0 = 0. 
which  is  i'.ot  very  accurate,  the  contribution  of  this  part  of  the  bod>  to  the  overall  pies- 
sure  drag  is  small  because  both/'  and  0 are  .small  near  the  point  of  maximum  thickiievs. 


l lie  pressure  on  the  afterbi  dy  from  the  position  of  maximum  thickness  to  tiie  ba.se  is 
ealeulaled  from  the  {’landtl-Mever  expansion 

i5i 

JO  y.M’'"  I 

( ate  must  be  exereiseil  to  ivstriet  the  negative  slope  on  the  aHerbrnK  to  values  levs  than 
about  S"'  i Ml  - ^ S" » or  else’  the  lltwv  will  separate  md  I <piaii*m  1 t w ill  m>t  Iv  valid,  * ' 

1 he  approximate  (uessute  distribution  given  by  1 1) nations  1 and  i5l  will  be  eompufed 
with  “exact'*  ttivts^id  mmteiieal  results  trotn  the  methrvj  ot  Solomon.  ei  al,'  ‘ 


SMh  I rivllon  Drag 

the  biiundary  layer  will  gHmeratly  K*  Uirbulent  4>ver  about  'Hi  jvrvent  of  the  proteetil' 
kHly  lor  kityte  caliber  ptoiecllles.  Sniev  the  lamimif  Mow  regum  i>  usuiilly  U>s  than  tU 
ivrevni  »H  the  total  surtace  area,  ti  will  K*  assmiievl  the  entire  iHunulitry  layer  is  uirhuleni.. 

( ruler  this  assninption  the  i«*tal  or  nieaii  skin  irietum  ovtlieu«l>  aevortlmg  to  Van 
Urtesi*^  may  Iv  obiai  . .)  from 

t/  7‘  » ' dstn  ' ^ sin  ' t .l 

M(\  t'  ' . 

wlietv  ■ 

21'  is 

£ ’ •«  - ..a:  - , it  W -■  . -w..  .*Y 

I//*  ♦ 4H’ * i/i''  r 4,1*  I*  ■■ 
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Uk'  variublc  n ol  I nualiun  (<>)  js  Uk‘  povvor  in  the  power  visei>sily  law 


ir^va 


ami  n tor  au  is  0.76.  hpialion  ((i.t  assumes  a tully  ik'velo|)ei|  turhuienl  tHnmUury  layer  with 
zero  pressure  gradient  ami  hanUtl  numlwr  ctpial  to  one. 

In  r*nlct  tr>  solve  I quatiun  t6i  tor  the  mean  skin  Inetion  coellieient  . one  must  have 
vaiiu's  for  / h‘?  I\  , A!  \ . and  . The  freestream  Key  nokls  immtvi  is  sirr.plv 

P r f 

»k  ♦ 


lo  relate  / / ^ to  the  ‘tvesireatn  Maeh  numK*r.  assume  the  wall  ;s  adial'alte.  Ik'lininga 

iMfhuleni  fi^wvty  favtorA*^  tty 


A',  •(k 

Ml  l»A/i 


it  is  well  knowtt  that  the  reetwery  laeitu  varies  a*  the  etihe  6H»t  ♦»!  the  l^randtl uunilvr  Ittr 
liithMUnit  flow  S4i  lh,i* 


Usvail  that  Van  Ihriest  V meth  yl  avamus  a rtaiutll  nuinlvf  »n  wmiy . If  this  \al«u*  were  usedv 
ll»s*ii  /*r  svouhl  also  (v  mi(i>  MowsWs'r.  Iln*  .ivtnal  sahu-  of  /’  ^ fl.7  f *o  tlt-il  the  ptevions 
ass«‘tiplt*nt  stf  haiidti  liniislvr  otie  van  ke  iom|teftsits^d  lor  sttmewhat  by  ths^  alnne  rsvtnvty 
faelot  whisli  foi  /^,  = 0.?^  vMHtUI  be  lltiis  I spiaiimi  l‘f|  itceomes 

vUi 


Then,  tor  a given  set  of  rreeslrcam  comlitirns  (M ^ ) one  cun  combine  ITiiuitions 
(-S)  uiul  til)  with  lb)  to  solve  tor  C.  . Tlic  eniuition  must  be  solved  numerically  however. 

’ oo 

since  (’,  cannot  be  solved  for  explicitly.  A procedure  adaptable  to  equations  of  this  type 
is  the  well-known  Newton-Kaphson  method. 

Once  the  mean  skin  friction  coefficient  lias  been  determined  for  a given  set  of  freestreani 
conditions,  the  viscous  drag  coefticient  is  .simply 


(12) 


1 he  wetted  area,  .V^, , is  the  total  surface  area  of  the  body  (excluding  the  base)  which  can 
be  integrated  numerically  given  a set  of  body  coordinates. 


base  Drag 

Much  theoretical  work  has  been  performed  to  predict  base  pre.ssure.''  There  is  still  no 
satisfactory  theory  available,  however,  and  the  standard  practice  has  been  to  use  empirical 
metlKHl.s,  Thi.s  is  the  approach  taken  here.  Tigiim  2 is  a mean  curve  of  experimental  data 
taken  from  Uefeamce  i 2.  Tliis  data  assumes  a long  cylimlrieal  afterbixly  with  a fully 
developed  turbulent  boundary  layer  ahead  of  the  base.  There  could  be  deviations  from  this 
cum'  due  to  low  body  fineness  ratio,  boattails.  angle  of  attack.  Keynokis  number,  and 
surface  temperaliire.  Tacli  of  these  effects  are  discu.sseil  in  Kel'erenee  12. 

I'he  ba.se  pressure  is  significantly  altered  by  the  presence  of  a hoaliail  .so  that  this 
change  must  he  accounted  for.  Probably  tlio  most  simple  melhoil  to  do  this  is  an  empirical 
cijualion  given  by  .Stoncy.'- 


and  thus 


(15) 


T.quiition  ( 1.^)  can  be  used  throughout  the  entire  Mach  number  range  where  ^ is  the  baa* 
passuiv  given  hy  tlie  curve  ofl'igure  2. 


OPTIMIZATION  METHOD 

The  optimal  sliapc  for  a fixed  length,  diameter,  and  Mach  number  is  determined  by 
computing  the  minimum-drag  body  fora  sequence  of  values  of  the  lorebody  length 
(refer  to  Figure  I).  I hen.  the  body  sitape  corresponding  to  the  value  of  which  gives 
minimum  total  drag  is  the  optimal  shape.  Hence,  two  distinct  optimization  problems  must 
be  considered:  the  i'orebody  and  the  afterbody  problems. 


Forebode 

A very  etficient  method  is  described  below  for  computing  the  optimal  forebody  shape, 
t he  algorithm  is  based  on  Lagrange  duality  theory  for  convex  control  problems  developed 
by  Hager  and  Mitler.  By  solving  an  initial  value  problem,  the  optimal  dual  multiplier  is 
computed.  The  optimal  shape  is  determined  using  a minimum  principle. 

,A  common  procediuv  for  solving  constrained  variational  probluns  is  the  following: 

An  iterative  metlunl  is  used  to  approx intate  the  solution! st  of  die  I’uler  ineiiiiality  associated 
with  the  vaiiational  problem.  I'ypical  algorithms  involve  Newton's  method  or  a shooting 
ineltuHl  to  solve  a two-point  boundary  value  problem.  Those  methods,  however,  am 
subject  to  ttumerical  instabilities  and  the  initial  iteration  imist  be  chosen  close  to  the  optimal 
siiajK  lor  ctmvergence. 

The  present  algorithm  is  far  Iviler  than  the  iterative  approaches  sinev  it  is  essentially 
direct  no  iterations  or  initial  guvsescs  are  involved  and  the  scheme  is  completely  stable. 

The  computational  procedure  is  described  Iviow. 


Kite  lombody  optimi/ation  problem  is 


minimiA' 

f' 

('  ftvtw'l  +f*tv)>  + 

.f‘(v)' 

1 

1 +f  1 

,vl^ 

• 

subjecUo  rtOi  « o.  rH:,»  = H 


r'tvl  > 0 (brail  ,ye|0,«^| 

where  tt  and  V,  an*  the  given  radius  ami  length  t)f  the  IbrelHHly.  Htc  constraint  r' t v)  0 
is  imp**setl  since  ft  v » Is  nondecreashig  on  llic  fomlHHly.  Aliluntgli  we  constrain  f tOI  » 0.  a 
bluni  n»tH*  shajv  ean  iv  aeliK*vs*d  as  iIh*  limit  of  a setpiertee  of  sluipv*s  sattsiying  ftOl  •»  0. 

I he  lihml  nose  slia|K*  emdti  alM^  be  treated  by  adding  an  additional  drag  lertti  to  the  eximmand 
above  iliai  eorres|H'nds  l<»  t|»e  nuliusof  the  nose  and  mnuKing  tlie  constraint  fit))  ^ Ot 
however  the  aiiaiysts  is  easier  using  the  tbnmilatioii  of  Lqiiation  ti  5)  above. 


H 


Make  a cliangc  of  variables  and  express  the  variational  problem  in  terms  of  the  inde- 
pendent variable,  p = r , and  the  dependent  variable,  j?(p)  = 1/r'.  Omitting  the  algebra, 
Hquation  ( l.S)  is  equivalent  to 


minimize 

1 

c>  VSpFTTt  /” 

5(P)2  + 1 

(16) 


subjeet  to 


g(p  ) 
s/p 


(ip  = 


gip)  ^0  for  p > 0 


lf,e(p)  solves  ‘iquation  (16),  then  the  corresponding  shape  r(.v)  is  found  by  integrating  the 
ollooing  initial  value  problem 


r'(,v)  = 

;•»())  = r„  (17) 


where  is  the  radius  of  bluntness  for  the  nose. 


As  shown  in  Kefeanta'  16.  twiulifferenl  tyi*es  of  lorebody  shapes  cun  occur  depend* 
ing  ot)  the  ratios(‘,,^^/C,  am'  iy</,.  The  optimal  shape  consists  of  a blunt  nose  followed 
by  a sinooth  monotoi'e  profile  whenever 


/: 


0 ■'*  ' 


(IH) 


If  this  inequality  is  violated,  the  optimal  sh:qv  consists  i f a spike  followed  by  a cone  with 
slo|H*  given  by  , 


.'(V.  . • i]' 


t or  realistic  tlighi  condiiions,  the  skin  Irictio  t crx’fficient.  (y , is  much  smaller  than  the 
stagnation  prv'ssuiv  c.x'fncient,  : hence  inetpiitity  ( IK)  is  satisla  d and  the  nose  is  blunt. 

The  algorithm  for  computing  the  blunt  nm>e  solution  oi  I qu.don  (16)  is  now  desertlvd. 
IKTine  the  following  futtciions 


*) 


1 


aiui 


//(ij.  r)  ~ //!> ) + - 


( l‘J) 


C0» 


riK'  tuiK'tion.  //.  is  ivlatcd  to  the  intcgraiul  of  ihi-  l.agraiittiaii  of  lii|uation  ( |t>)  j^ivcii  hy 


r V(P)'|  f'^‘  / 

vP  J \ 


//|V(P).  cJp 


By  tiu'  miniimim  priiuiplc.''’  if.c’^tp)  is  an  optimal  solution  to  1 iiuation  i IP).  iIkto  exists  a 
scalar  X such  that 


//{je*(pK\/p/X)  - niin{//(p'.vP/X):  a > Oi 


for  almost  every  0 p < A'* , lienee,  after  eoniputinjt  the  optimal  ilual  multiplier 
X..e*  (p>  is  the  value  of/»  > 0 that  minimizes //(^>.>/p7X). 

Let  tH:  I he  the  value  of  *•  P 0 that  mininti/es  //(jl*,;  ).  AnalyzitiM  the  structure  (ft  the 
fimetion  Aljel.  it  ean  Ik  shown  that  there  exists  a critical  value. ; such  that  (/t;)  = 0 
lor : ^ and  (i{;  l > 0 for  .*  > ;,j : this  critical  value  is  given  by 


:» * t/t'tv'ir"  T))’' 

wheix*  0 is  the  uni(|ue  positive  rmn  of  the  etpiation 


Cl) 


122) 


Also,  it  can  be  shown  that  tlie  slope  of  the  nose  at  \ « 0 is  l/v'SF'^n”  and  the  slope  upproaches 
45“as(‘  /<),  -0. 

Let  u)U  ve,  siitisfy  the  dil'I'eK'nUai  ci|ualions 


te', (:)  e (#■(:)  n ^t;)  «*  2yti^ 
te I fr»  I ^ t>  u* j I- 0 J “ ^Uo L',  \ld, 


C3f 


Notice  il).)t  n u^(;„)a«ul  as;  “*  *. <;t; I approaches  the  left  side  of  inctjUiiliiy  ( IK). 
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Tims  by  inoiiuality  ( 18),  u'J(r)  > >»',(r)  for ::  sufficiently  large,  there  exists  5 such  that 
u'l  (5)  = w'2^-)-  It  can  be  shown  that  the  optimal  dual  multiplier  is  X = f|/f  and  the  nose 
bhmtness  is  - X:„.  tlenee.  the  solution  to  I quation  ( Ibi  is  g*(r)  = and  the 

optimal  shape  is  given  by  I quation  (17). 


Afterbody 


In  addition  to  pressure  drag  and  skin  friction  drag,  the  afterbody  also  has  base  drag. 

I quation  (,S)  must  be  integrated  to  obtain  the  pressure  drag,  and  the  base  drag  is  given  by 
l-()uation  (Id).  The  afterlwdy  drag  eoelTieieiU  to  be  minimized  thus  becomes 

^/)  I ' (~  l]  r(v)r/(.v)  + (■  r(.v)v'T+</7  . 


(\,  fit!)' 
1)1  \ 

ikn 


subject  to 


r(.vt  ~ </(v) 


(/‘tv)  " tMv) 


, TfH(V)  /»(v).W(/>l* 

'"'’‘iTsv  'c^r 


The  initial  coiulitiuns  are  at  the  txgimting  of  the  aljcrbody  where 


(:4) 


CS) 


rt^^) «(/,/' 


i/tt(,)‘'  0 

/HV,)  - /». 


t:<>) 
(.'7 1 
(281 


I he  \ iiialde  ;i(  \ ) »/*(  v ) ^ / ’’( v)  is  t»*  Iv  elu»><*n  to  minum/e  tlte  drag  cwificieul : I v|uaiion 
(-7)  implies  the  IhhI>  slope  at  the  iH'gtnniiig  ol  the  aflerN*d>  is  zero.  \ovv  the  slo|v  at  the 
end  of  the  fim?t'"d>  is  generally  non/eto.  therel«»re.  there  can  Iv  a slo|v  dlseontutuity  at 
V t,  I quation  ( 'X»  gives  the  Ivginning  pressure  as U)  which  is  eottsisteni  with 
nunlified  Mew  tonian  ihemv , I Iv  Mach  number  on  the  surlace  is  related  to  the  surfasv 
prvssuiv  b>  the  isc‘niro|iie  relaiion 


« \ 


(.'*») 


M 


whore  the  stagnation  pressure  is  given  by 


/»o 


(30) 


In  i:quation  (35)  the  inite|)eiulent  variable  for//(.v)  was  changed  from  0 in  liquation 
(5)  to  V by  use  of  the  geometric  relation 

do  ti  ( V) 

— = (3 

l+</(.v)^ 

As  used  m control  theory  terminology,  the  fimctk>n  ii  above  is  the  control  variable  and 
the  functions n q.  and/>  are  the  state  variables. 

I he  I tiler  approsimation  is  used  tt>  replace  the  integral  and  derivatives  in  lii|uation 
(24)  by  the  lollowing  discixii/athm 


mmimi/e 


subject  to 


1 

r > .\ 

V 

\ 

I 

1 

1 

" 1 ^ 

, V 

Ml 

til " ) A\i 

HI  \ 

1 

tt  I) 

J 

> 

*1 

«!.  1 ** 

+ A v /•'(  1 " . fi"  ) 

r**i 

V 

« a 

I'i 

m 

(/" 

► 

L//'J 

i4/2, 

. «1 

* 0. 

(32) 


(33) 


where  |r«.  t/",  />"  | aiv  appruNimalioits  (o  K ).  /»( '^  M . q.  /»» deiipies  the 
tniegnuid  in  Ikiualioii  (24).  and  Ihe  wclor  function  Fir,  q,  /».  n)  denotes  the  ri|dti  side 
of  the  dirfea*ntial  equations  in  I quution  ( 25). 


Fite  Mthilion  to  I quation  (32)  was  cont|>uled  by  (he  itH'thod  of  sbvpest  descent. 


IX*  line 


//( I " . . r ) Lu  " ) ♦ / ( t “ . II"  )♦  X" 


(.34) 


and  let  V < denote  iisf  ol  dual  ntulliplk*r  vectors  gnen  by 

‘l  ■ 

y ' s r ♦ Av  //( « '*  til"  I.  fi".  V’ ) 

hi  « 


(.351 
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where  \"(ii")  is  the  solution  of  tlie  ilitferenec  I'quation  (32)  for  given  {ii"  }.  Starting  with 
an  initial  guess  to  tlie  optimal  control,  tire  control  is  incrementeil  by  the  rule 

=■  *Cm,  - (361 

air 

I or  the  argument  of//  in  btiiiution  (36),.e"  is  determined  by  the  difference  liquation  (32) 
and  X"  is  determined  by  l-.quation  (35)  for  given  The  parameter  f > 0 is  chosen 
small  enough  to  reduce  the  drag  for  the  new  disease  control  approximation.  Hor  the 
bi>dies  studied  in  this  report,  c % 40  gave  the  most  rapid  convergence  in  the  discrete  control. 


KliSULTS  AND  DISCUSSION 


COMPARISON  OF  PRFSSUKF  PKFDICTION  METHODS 

To  test  the  accuracy  4 the  methods  used  to  calculate  the  pressure  drag.  I (|ualiuns  (2) 
and  (.5)  are  compaa'd  whb  numerical  a‘sulls  from  the  method  of  Solomon,  el  al.' 

I igua's  3. 4.  and  5 compt  v the  passure  distributions  over  the  original  nose  boattail  body 
given  by  the  equation 

ra|l,^tlv  v^)'*  y/lT.  0 Sss  V 1.75  (37) 

for  .1/^  s 2.  3i  and  5.  a speclively . I he  toreNnly  pre«kure  was  computed  by  two  methods: 
ta)  F(|uation  (2).  and  *b)  the  mtHlilled  NewtoniamBitsemann  passua  distribution  which 
includes  tlw  eurvatua  effect.*  These  three  Hgua*s  show  that  I quation  (2)  padiets  the 
foabody  pressua  aaamably  well,  and  Fquation  (5)  gives  panamas  aostmably  elose  to 
the  ‘V\act“  solution  oh  the  afterbtHly.  .Actually  tia  calculated  drag  is  quite  goitd  Kaause 
the  passtire  is  trm  hiph  near  the  maximum  thickness  and  imt  low  near  the  base.  These 
two  effects  intaxiuev  partially  comiHmsating  errots.  In  all  three  figuas.  the  cunatua  effect 
drives  the  passua  t)H>  low.  hence,  this  method  was  not  used  in  the  optimi/ation  pixKVSs. 


OPTIMAL  SHAPES 

l or  fixed  Vidues  of  t U and  3/^ . minimum  drag  Nnly  shapes  wea  calculated  tor  foa- 
IwHiy  lengths  of  A « 0.5. 0.b.  0.7. 0.H.  and  O.M.  | l;en  the  itptimum  value  of  V (for 
a giveit  T:t/  and  3/^  f is  tlw  value  eorrexponding  to  the  minimtimC /) . Eigures  6. 7,  and  b 
sliow  the  efieci  of  V on  ( fw  values  of  4. 5,  and  C».  aspeetively . Each  figua 


uivvs  till.'  resiilt.s  t'or.V/.,  = 2,  3.  4 and  5.-;  Noif  tluil  each  curve  indicates  that  the  '^altie  of 
t,/\’  corresponding  to  nuidnuim  (.'0  is  not  sliarply  del'ined.  In  fact,  Figtiie  9 show.s  that  for 
Kjd  - 5,  which  gives  minimum  varies  almost  linearly  with  , and.  that  AV^/6  ^ ± 0.05 
only  increases  the  minimum  Cj^  by  about  ] percent.  Also.  Figure  10  illustrates  the  effect  of 
on  ( for  = 3 and  C/t7  4,  5.  0,  and  7,. 'A  cross  plot  of  this  figure  is  shown  in 
l-‘igiire  I I which  depicts  the  optimal  as  vaiyirig  almost linearly  with  i/d  and  a range  of 

Af  ./t:  ==  1 0.05  onlv  increases  C/)  by  about  1 percent. 

^ m i n 

File  optima!  shape  for  = 3 and  tVd  = 5 is  illustrated  in  Figure  12.  This  shape  is 
: close  to  that  found  by  Moore'  using  scmiempirical  techniiiues.  The  forebody  has  a value 
of  0.7  and  the  ba^c  diameter  ratio  hd^Jd^  = 0.7.  From  consideraiions'of  the  internal , 
ballistics  of  projectiles,  the  value  of  ty/t!  is  limiTcd  to  about  0.05  to  0.7  unless  forward- 
mounted  sabots  arc  used  for  launching' the  projectile.  However,  as  shovyn  in  Figures  9’arid 
1 1,  even  though  the  point  (yf  ina.xinnim  diameter  inay  not  be  at  the  optinium  location, 
only,  a .snudFdrag  [k'halry  is'paid  for  a ± 5-percent  variation  in  iji  from  the  optimiim. 

. fjgure.  13  shows  the  variation  of  f,,  with  for  t/d  = 4,  5,  and  6.  The '’iirve  desig- 
nated C is  the  drag  coefficient  lor  the  optimal  shape  at  each  A/  and  ijd.  Hence  the 
'mm  . 

value  for  is  different  for  each  A/^  and  i/d.  F' .'  dashed  curves  give  the  variation  of 

w'ith  A/^'  for  a fixed  value  vXijIi  = 0.7.  This  figure  illustrates  that  = 0.7  gives 
nearly  inininium  drag  for  i/d  = 5 and  6 over  tlic  range  2 < A/„  <5. 


CONCLUSIONS  AND  RECOMMENDATIONS 

'.  A numerical  method  has  been  developed  for  caieulating  the  optimal  projectile  shape 
which  produces  minimum  total  drag. 

2.  The  optimal  shape  changes  withA/d  and  Mach  number.  I'or  i/d  = 5 a nose  length  of 
i^/i  - 0.7  gives  nearly  the  uptiimiin  shape  over  the  raiige  of  2 < A/^  < 5. 

3.  A variation  in  no.se  length  of  MjH  J 0.05  results  in  only  about  u l-|H‘rccnt  inewase 
in  total  drag. 

4.  The  optimum  .shapes  have  a base  diameter  of  about  0.7  limes  the  maximum  diameter. 


5.  It  is  recotumeiuleil  mat  a similar  optimization  Ix'  performed  for  subsonic  and  tran.sonic 
Macii  miinbers. 


(1,  While  it  is  holieved  the  optimum  body  shape  tierived  from  the  approximate  theory  is 
reasonably  aeeurate.  shotiM  an  exact  theory  he  used  for  the  pressure  prvdietion  instead 
of  the  approximate  tlieory,  the  absolute  value  of  the  drag  etx'fndent  may  Iw  in  sub- 
stantial error,  particularly  at  the  'ow  supersonic  Mach  numbers.  For  this  reason,  it  Ls 
believed  Ibe  expenditure  involved  in  using  an  exact  theory  such  as  that  of  Solomon*-' 
to  repeat  the  preseu*  work  is  jiistitied.  Tilts  effort  will  he  started  in  l-'Y  77  at 
N.SW(7I)|.. 
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APPENDIX  A 


GLOSSARY 


lunctions  (icfincd  by  Lquation  (23) 
distuncv  along  body  axis  of  symmetry 
parameter  used  in  Equation  (20) 
critical  value  of  : given  by  Equation  (21) 
positive  root  of  Equation  (22) 
ratio  of  specific  heats  (1.4  for  air) 
angle  along  body  surface  d = tan  - 1 iJrldx) 
dual  multiplier 

cwfficient  of  absolute  viscosity 

air  density  also  p = in  forebody  optimization  problem 


GLOSSARY 


A 

B 


Co 

(\ 


AB 


c. 


c. 


pba 

pba 


^B 

»• 

a 

it 

H 


» 

P 

/*» 

A 


constants  donned  after  Kquation  (6) 

drag  coefficient 
afterbody  drag  coefficient 
forebody  drag  coefficient 
mean  skin  friction  coefficient 
pressure  coefficient,  c'^,  = p ~ 
base  pressure  coefficient 

base  pressure  coefficient  on  a configuration  with  a cylindrical  afterbody  and 
turbulent  tlow  at  the  base 
stagnation  pressure  coefficient 
base  diameter 

reference  diameter  (maximum  cross  section  of  body) 
l/r* 

optimal  solution 

value  of  g which  minimises //fg.;) 
function  defined  by  Lquation  ti^) 

function  defined  by  liquation  (20).  abo  a different  function  by  liquation  (M) 

length  of  configuration 

forebody  length 

local  Mach  number 

frecstream  Mach  number 

variable  defined  by  I qtiation  (7) 

presstire 

slagnat  km  pressure 
^ndtl  number 
freesiream  pressure 
maximum  radius  of  body  « </,  /2 
freestream  Reynolds  numK*r 
turtiUlent  recovery  factor 
radius  along  body 

blunted  nose  radius  at  lip  of  forebody 

body  slope. 

reference  area  * » 

weItvM  surface  area  of  botly 

wall  temperature 

freestream  tenvperattire 

freesiream  velocity 

A*l 
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